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ABSTRACT 


Submersibles are limited in depth by the strength of the pressure 
hull. The pressure hull, in turn, is limited in strength by the properties 
of the material used in fabrication, construction techniques, and by the 
weight of material used in construction. Once the ballast weight has been 
taken up in hull weight, a limit is reached in depth. To go deeper with the 
Same _ payload, the hull must be enlarged or size may be held constant, the 
hull strengthened and flotation material attached to support the extra 
Weight. Either of these methods, enlarging the hull or adding flotation, 
increases size greatly as depth is increased. 


This thesis is a study of total vehicle size versus depth for 


weight-limited deep submersibles with and without buoyancy material 
attached. 


Thesis Supervisor: Philip Mandel 


Title: Professor of Naval Architecture 
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Ie DISCUSSION 


Ae Volume-Limited Submarines 

1. Conventional submarines have operated routinely 
in the upper layers of the ocean for the better part 
of a century. These submarines, designed primarily 
as military vehicles, are large, providing the neces- 
sary internal volume required by men and equipment. 
Weight of hull and equipment is less than the sea 
water displaced, thus it is necessary that these 
vessels carry ballast to attain neutral buoyancy. 

2. To design a submarine for maximum depth, it is 
necessary that the hull weight be increased. This 
can be accomplished with a volume-limited submarine 
with no increase in size by trading ballast weight 
for hull weight. The depth to which a volume-limited 
Submarine can be designed with no increase in size 
is then, a function of the ballast 1t carried. Beyond 
that depth the weight of hull and equipment exceeds 
the weight of sea water displaced, and the vehicle 
becomes weight-lLimited requiring additional buoyant 


support to maintain neutral buoyancy. 





Be Deep Submersibles 

lL. Oceanographers require small laboratories to carry 
them to the depths of th sea. The desigmer's problem 
in building submarines to satisfy the needs of ocean- 
ography is to build vehicles capable of carrying a 
Scientifically useful payload to a scientifically 
interesting depth, usually the maximum depth possible. 
submarines of a given size, Limited in depth by 
hull weight and hull materials, have several avenues of 
approach by which their depth may be increased. First, 
vehicles may be made Larger on payload fixed. 
Improved pressure hull materials may be used or flota- 
tion material may be attached to the vehicle to provide 
the lift necessary to restore neutral buoyancy. The 
first and last of these methods result in large craft which 
are less maneuverable, more complicated, and more expensive 
to build. This thesis examines and compares the penalties 
paid in achieving maximum depth by the methods just des- 


cribed utilizing various pressure hull materials. 








il. PROCEDURE 


A. General Approach 

i. A number of assumptions are made concerning the 
family of vehicles to be considered, the structural 
materials available, the desired payload, the pro- 
pulsion characteristics and the other features which 
Will affect the buoyancy-welght balance. The 
assumptions are not necessarily optimum for the range 
of vehicles considered, but they represent an average 
of the features of TRIESTE II, ALUMINAUT, and ALVIN 
as they apply to the family of vehicles considered. 

2. <A computer program was developed to optimize hull 
weight for a given depth. The program is described 
in Appendix A. Using the computed results of hull 
weight versus depth as a base, the weight of payload, 
batteries, dropable ballast, main propulsion units, 
and "all other" variables was added. This resulted 
in a basic vehicle capable of proceeding to a given 
depth without the need of flotation material. At 
Lesser depths additional ballast is required and at 


greater depths added lift is necessary. 








3- Vehicles of various size were used in the computation 
of hull weight versus depth. Size was varied in two 
ways: by holding outside diameter constant while increas- 
ing length and by holding the length to diameter ratio 
constant. The results determine changes in size with an 
increase in depth for the two methods of increasing dis- 
placement. 

4. Finally, flotation is considered as a means of hull 
Support at increased depth. The addition of flotation 
necessitates an addition of the variables mentioned 
above, batteries, ballast, etc. This reduces the effect 
of the flotation material. Curves were constructed for 
Syntactic foam, probably the best flotation material now 
available, and for materials of lower densities which 
might be developed in the not-too-distant future. These 
curves take into account the weight of additional equip- 
ment to arrive at a reduced flotation effectiveness. 
Entering the curves with the additional lift required by 
a basic vehicle at a given depth, one determines the 
volume of flotation material required, weight of flotation, 
and the weight of the added hardware necessitated by the 
additional displacement. The result of primary interest 
is the displacement of flotation to be added to the dis- 
placement of the basic vehicle. Total displacement is then 
plotted against depth for comparison of size with similar 


curves for the larger pressure hulls considered. 


Jee 








B. Assumptions 


A family of cylindrical vehicles with hemispherical 
end caps is considered. 

A basic vehicle with outside diameter of seven feet 
and 1i/D ratio of 2.0, where Lis the Leneuh smeindin 
the hemispherical ends, is considered. Hulls with 
eight-foot and 10-foot diameters, also with L/D ratios 
of 2.0, are considered for comparison as are vehicles 
of seven-foot diameter with L/D ratios of 2.5 and 3.0. 
Thus, size is varied in two ways; by holding L/D constant 
anid oy eoolise eons Cant. 


Four materials are considered for construction as 


follows : 
Material Yield Strength Modulus of Elasticity Density 
Steel SO,000 pace 30 x 10° Tececee 4.90 lp. /ft.> 
Steel 508000 Wee. i. penx 10° eee 4.90 lb. /ft.> 
Mi taniun aC OO Das anes ee 10° pos sale 281 Lb. /ft.? 
Aluminum GO COOmpes: is 10 x 10° secs 168 1b. /£t.> 
Glass was also included for comparison in the computation 
of hull weights, using the same theory applicable to metals. 
This theory may, but probably does not give realistic values 
for the weight of glass hulls. The properties used for glass 
auerao HOLLows : 
Glass . 300,000 p.sei. 14 x 10° foe Sic 184 lb./ft.> 








The values of hull weight for glass are tabulated 
in Tables I - ITT but are not plotted. 

4. Power requirements are based on a sustained speed 
or five knots for a period of ten hours. The auxiliary 
load is assumed to be 1O kilowatts. 

5- Main propulsion power is assumed proportional to the 
product of displacement to the two-thirds power and the 
speed in knots cubed. Since the speed is fixed at five 
knots, power is proportional to displacement to the 
~wwo-thirds power. 

6. Silver-zine batteries with the following properties 


are assumed for the power source. 


Capacity 6.5 kilowatt hours/cubic foot 
23.0 pounds/kilowatt oe 
Te Payload including complement and instruments is held 
constant at 2,000 lbs. for all vehicles. 
8. Minimum volume requirement for payload and equipment 
was assumed to be 366 cu.ft. which is the internal 


volume of the seven-foot diameter, L/D = 2 pressure hull. 





(1) Keays, K., “Parametric Study of Two Man Search and Work Submarines," 
General Dynamics/Electric Boat contract study, page 24. 





9. Auxiliary buoyancy is assumed initially to be 
syntactic foam with a density of 44 lbs. per cubic 
foot. Buoyant materials with Lesser densities of 
0.50 (32.1). per cu.ft.) and 0.30 (21.33 per cu.ft.) 
are Studied for comparison. 

LO. The remaining weights which do not fit into the 

categories of hull, batteries, or payload vary with 
a number of parameters. Some vary with displacement, 
Some with depth, some with speed and some with external 
or internal surface area. It is not the purpose of 
this paper to analyze each of these welght groups, but 
estimates must be made in order to concentrate on the 
buoyancy requirements. It is assumed, therefore, that 
the following weight categories in addition to hull, 
battery and payload will be considered. 

Main Propulsion units - function of displacement - 

3% of displacement 


Flotation material - function of buoyancy requirements 
and depth 


Ballast ~- function of total displacement - 8% of displacement 


All Other - for the purpose of this thesis a ns 
of displacement - 8% of displacement 2), 





(2) Wenk, E., Jr., DeHart, R. C., Mandel, P., and Kissinger, R., Jr., 
"An Oceanographic Research Submarine of Aluminum for Operation to 
15,000 Feet", Transactions R.I.N.A., Vol. 102, No. 4. Oct. 1960. 
Based on estimates made for ALUMINAUT, page 15. 
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C. Basic Vehicle Weights 

L. A computer program was developed based on the class 
notes from course 13.15, M. =. T., by J. Heme 22 4a-. 
The program optimizes frame size and frame spacing to 
arrive at a Least weight solution for a hull of a 
given size, constructed of a given material. The com- 
puter program is described in Appendix A. The results 
are tabulated in Tables I - V. 

2. Figure VI is a graph of battery weight versus total 
displacement. For vehicles above 30,000 lbs., battery 
welght can be treated as a constant plus a percentage 
of the displacement over 30,000 lbs. with an error less 
than one-tenth of one percent of total hull displacement 
at all times. For a vehicle to which flotation material 
is to be added, the additional battery weight equals 
1-4 per cent of the additional displacement. 

3. Eallast, propulsion units, and "all other" also are 
functions of displacement and, like batteries, these 
weights may be treated as a constant which is part of 
the basic vehicle, plus a weight equal to a percentage 
of the displacement of the buoyant material added. Added 
weight due to an increase in displacement amounts to 


20.4% of the added displacement. 





4. Figures II - V are plots of weight versus depth and 
take into account hull weight, batteries, ballast, 
propulsion units, "all other", and payload. Weights 
plotted in Figures IIT - V are exclusive of any flotation 
material. Displacement is also plotted in FiguresII - V 
for the differently configured vehicles. Comparison of 
displacement to weight yields the ballast or the added 


Lift required for neutral buoyancy. 


D. Buoyancy 

1. The buoyancy of a vehicle is equal to the weight of 
salt water displaced. Flotation material must be added 
to a vehicle which weighs more than it displaces; ballast 
must be loaded aboard a vehicle which weighs less than 
its displacement. 

2. The addition of flotation material necessitates an 
addition in the requirements for batteries, dropable 
ballast, main propulsion units, and "all other". Having 
linearized battery weight as a function of added dis- 
placement, all of the above items may be treated as a 
Single weight addition proportional to the flotation 
material displacement. 

3. The lift produced by flotation material is also a 
percentage of the flotation material displacement. The 
percentage may be treated as the effectiveness of the 


material. Syntactic foam, at 44 lbs. per cubic foot 
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is 31.3% effective. <A material with a density 
with respect to sea water of 0.5 is 50% effective 
and one with a density of 0.3 is 70% effective. 

4. The effectiveness of the flotation material is 
reduced by the additional weight described in 
paragraph 2 of this section. The reduced effective- 
ness is the flotation material effectiveness less 
the percentage of added weight. For the assumptions 
made in this study, the reduced effectiveness of 
syntactic foam is 10.9 (31.3 - 20.4);of the 0.5 
density material it is 29.6; and of the 0.3 
density material it is 49.6. Since this reduced 
effectiveness is the true indication ox Littecentri. 
buted to a basic vehicle, it can be seen that | 
reducing the density from 0.687 to 0.5 increases the 
effective lift almost threefold. Flotation material 
displacement is reduced to one-third of the initial 
requirement and added batteries, dropable ballast, 
etc., is only one-third as much. 

5. Figures VII - IX are graphs of flotation material 
characteristics and weight additions. Entering along 
the ordinate of the graph with the lift required by a 
basic vehicle at a given depth, the curves provide the 


following information: 


2 One 








Curve 1 


Curve 2 


Curve 3 


Curve 4 


Curve 5 


Lift applied to the basic vehicle 


Lift applied to additional batteries, 
ballast, etc. 


- Total lift provided by the flotation 
material 


- Dry weight of flotation material 
to be added 


- Salt water displacement of flotation 
material to be added 


These curves apply to any of the vehicles in the family 


of vehicles considered, provided that flotation material is 


to be added to provide additional lift. 


oe 





TIT. PRESENTATION OF RESULTS 


Pressure Hull Sizes 


Figure I 
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Displacement of Hulls 


Hull Diameter Hnre@ ats L/D Displacement 
(Lbs. Salt Water) 
ie ape 230 28 ,650 
Be aS aS Bi) 300 
ieee 3.0 45,950 
SC LCRO 2.0 42,900 
Oe 2On0 Za 83,500 
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Table I 


Computed Cylinder Weights L/D = 1.5 


Collapse Depth 


Ee ere eS 
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Aluminum 60 





Glass 
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Table cs 


Computed Weight of 2 Hemispherical End Caps 


Collapse Depth 


| By OO. fice MOpe CO ance if O0Oe fcr 3) COG asic 
Steel, HY-80 | 
2,300 lbs. 7, 3, 00Mlpcr i OO Cnn ae 22,700 lbs. 
Pe] ano | some | anew | ssom 
7,200 21,200 34,700 66,100 
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Table III 
Computed Pressure Hull Weights L/D = 2.5 


Collapse Depth 
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Table IV 


Computed Pressure Hull Weights L/D = 2.0 


Collapse Depth 


> Be hOeraiue 10,200 .Gts LOO cs 34 000crt. 
Steel, HY-30 

fer 9,100 25,400 78 , 700 

17,600 49,000 79,200 51 LOO 
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Table V 


Computed Pressure Hull Weights L/D = 3.0 


Collapse Depth 
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Pl == a 
ee 
pel eae oem 


oan .um 






















Aluminum | 


pr] sco || sco | sac 





6 , 700 17,500 28 ,800 54,800 





ae 





Table VI 


Weights Added to Bare Hulls +o Arrive 


Diameter 
L/D 


Battery Wt. 
(Figure VI) 


Ballast, Etc. 
(19% x Displ. ) 


Payload 


Ocak 


L/D 
Battery Wt. 
Ballast, Etc. 


Payload 
Total 


at Basic Vehicle Weignts 


ie 
2.0 


3,300 


5,450 


OES 


LO, 150 


a 
2.0 
3,600 
ope 


2 , 000 


13,130 


(Figures II-V) 
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Figure XII 


sample Computer Solution 


Material - Aluminum 60 


Collapse Depth - 17,000 feet 


Hull Diameter - T7' 





One 





IV. CONCLUSIONS AND RECOMMENDATIONS 


A. Figure X displays the results required by section I of this thesis. 
Overall vehicle size is plotted against depth for hulls constructed of 
four different materials. The effect of increasing depth by changing 
hull size or by adding flotation material with pressure hull size con- 
stant is compared with the following results. Increasing hull size by 
varying vehicle length provides greater depth than building a larger 
vehicle with L/D constant for HY-80 and Titanium. For aluminum the 
opposite is true and for HY-150 both methods have equal effect. Therefore, 
the results with respect to a preferred method of increasing size are 
inconclusive, one method having nearly the same effect as the other. 
However, the addition of flotation material is to be preferred over 
changing vehicle size by either method for all but HY-80, assuming 
Syntactic foam is the material used. For depth increases greater than 
Emeguc 5,000 feet, a smaller vehicle results if flotation material is 
added. And, of course, the cost of a given volume of flotation material 
will be mich less than the cost of increasing hull size by an equal 
amount. 

Recommendation: Deep submergence pressure hulls should be 
built large enough to satisfy minimum volume requirements. If added 


lift is necessary, it should be provided by flotation material. 
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B. The effect of improving the characteristics of flotation materials 
is graphically demonstrated in Figure X. If a flotation material could 
be preduced with a density of 0.50, there would be no question about the 
preferred method of increasing depth. 

Recommendation: As flotation material of lower density is 
developed, its use should be considered in the design of deep submersibles, 


Costs will be cut and smaller, less complex deep submersibles will result. 


C. Although the results concerning a preferred method of changing vehicle 
size were inconclusive, it can be said that there is probably an optimum 
fee contipuration. With the computer program developed, one could azrive 
at hull proportions for a least weight vehicle with a given displacement. 
At first glance it might appear that this would always result in a sphere, 
but with conventional framing or sandwich hull construction, this probably 
will not be the case, The program can also be used to determine optimum 
frame configuration, and it appears to be a very promising tool in evalu- 
ating current design practices as they apply to vehicles designed for 
greater depths. 

Recommendation: The computer program described in Appendix A 
should be used to continue the study of hull configuration started in this 
paper and it could be used as well for a study of the parameters of sub- 


marine design and their interaction on one another. 


D. Figure XI is a graphic illustration of the results obtained by using 
the methods outlined in the notes from course 13.15 and from earlier cal- 


culations based on unstiffened cylinder construction. Some difference 


orale 





exists between the weight of unstiffened cylinders and the weight of 

the stiffened cylinder as predicted by the program in Appendix A. 

It is natural that a difference exist, but in some cases the unstiff- 
ened cylinder is limited to a depth 38% less than that of a stiffened 
Gyiinder of equal weight. This large difference is difficult to explain 
and it may be an indication that the notes on which the computer program 
is based may give incorrect collapse depths when considering deep sub- 
mersibles, or, that the computer program contains some flaw. However, 
the program was tested using as inputs data from relatively shallow- 
depth submersibles, and it duplicated almost exactly the calculations 
which resulted when those vessels were designed. 

Recommendation: The computer program, probably the most 
valuable by-product of this thesis, should first be checked once again 
to ensure that the formulae have been translated into computer language 
properly, and then the output should be compared with experimental results 
to determine: 


Ll. the validity of the theory with respect to 
conventional steel construction, 


2. the applicability to other metals, and 


3- the application of the theory to vessels 
with ce ratios approaching 0.1. 


Having checked the program and its general applicability, it will be 
possible to determine the necessity of modifying the basic theory for 
metals other than steel or, more probably, for deep submergence vehicles 


With very thick shells. 
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APPENDIX A 


Computer Design Program to Optimize Submarine Hull Weight, based on 
elass notes of Course 135.15, M. I. T., by deere, byvancee 2c. -. cu 


of Naval Architecture. 


ime background 

a. Failure of a stiffened shell occurs in one of three 
modes. The first occurs if the shell buckles between 
the frames while the frames retain their circularity. 
This is referred to as elastic or plastic instability. 
The second mode of failure involves ylelding of the 
shell. In this mode the shell remains circular but 
the yield strength of the material is exceeded. The 
third mode of failure occurs when the shell and frames 
fail together and this is known as general instability. 
Ideally, for optimizing strength with Least weight, a 
shell should fail in all three modes simultaneously. 
This program optimizes three parameters, shell thiclmess, 
frame area, and frame spacing in an attempt to find the 
least weight solution. 

be In such a solution it would be ideal to optimize the 
dimensions of the frame as well as the other parameters 
mentioned, but this would have introduced another degree 


of complexity to an already complex analysis. Instead, 


aa 





a frame typical of those presently used in submarine 
construction was selected as a parent and all frames 
are scaled from the parent member. The cross section 
of the parent member is shown in Figure A-1 below. 
All dimensions are with respect to a web 


thickness. 


Shell 





ii 
Je 6.8, =e 


c. Having selected a parent frame member, it is only 
necessary to settle on the length of cylinder between 
bulkheads and the mean cylinder diameter in order to 
apply the formulae presented in the set of notes refer- 


enced above. The general method of approach outlined 
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is that of selecting a set of scantlings and checking 
heir applicability. An understanding of the paper is 
necessary for a full understanding of the program and 
its capabilities. However, lack of such understanding 
will not prevent successful use of the program. 
d. The program is based on the solution of six equations. 


These are: 


L- Hoop stress formula which predicts a pressure 


at which shellwyield oceurs. 
Per = 2(+/D) Oy, 


This formula is used to determine shell thickness 


for a contemplated design. Shell thickness, 


a Pex ae 


B 
oy 


7G 


Experiments indicate that the predicted thickness 
may be reduced in some cases by as much as 20%. 
The computer program, however, completes two 
solutions, at predicted 't' and at 95% of 't' but 
a minor modification to the program allows 't' to 
pe reduced to any given percentage of the value 


predicted by the hoop stress formula, or it allows 


= om= 





1+. 7 
u 


to be reduced in increments unvil a least 
welgnt solution is found wnich satisfies all 
other criteria but ignores the hoop stress 
predietson. 

2. Formla 92a modified predicts plastic insta- 


bility failure. 





a 
D. v 1 
Por * a a en 
i. 70:85 - B L + O.LU8(X = ¢) 
L+ 6 


The crivical pressure it predicts mist be equal 
to or greater than collapse depth pressure. (see 
page 42 for a list of symbols and their meanings. ) 

3. Wenk's formula predicts collapse pressure based 
on the circumferential stress at the shell mid- 


thickness midway between frames. 


eu 1 
P oa aue x 


er D 
» | 2 a 


kh. Lunchick predicts "yield pressure” based on a 


three hinge failure mechanism. 


ie eeeee Pe - 


ct Vie coi 


This formila is normally within 3 - 5% of the collapse 


pressure predicted by Wenk. 
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o« Bryant has taken a collapse prediction by 
Kendrick and put 1t inve 2 eerie case es 
Tormila is based on the number of complete cir- 
cumferential waves formed by the plating and 
frames buckling as 4 Unite iigeis solved emer a. 
for 2, 3, or 4+ waves and the first critical mode 
of failure (i-e., the one for which the critical 


load is smallest) is the value sought. 


h 2 
> = ent . m x 8(n - 1) Fr 
er D, Le 





(no + a 1) (ne + oe “n : 
Formulae 2, 3, 4 and three variations of 5, predict 
critical pressures. It is necessary to solve these 
formulae and then select that answer which is lowest 
as the critical one in the case under consideration. 

6. ‘The last formula used predicts the stress in the 
frame flange due to the combined effects of total 
load including shear and bending due to out of 
roundness. This formula is based on the critical 
pressures predicted by Bryant's formula} thus, it 


is also solved three times in the program. 


L + 0.858 
ai B = 
alk E(c ) ie i \ 
Cmax ~ e ae (n” = 1) i(5 = pe 
A, doce 3 - ex 
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The highest values predicted by this formula 

mist not exceed the yleld strength of the material. 
Thus, the highest value is selected by the program 
as the critical value. 

e. The program solves for shell thickness using the hoop 
Stress formula and then computes six solutions for critical 
pressure. The lowest of the six is termed the critical 
pressure for solution of the problem. ‘Then three solutions 
of maximum stress are made, and the highest of these is 
Selected as the critical stress. At this point the program 
computes the amount by which critical pressure must be 
changed to agree with collapse pressure, and the amount by 
which the highest stress mist be altered so it does not 
exceed yleld strength of the material. If the program can 
be made to alter frame size and frame spacing until these 
values agree; critical pressure with collapse pressure, and 
maximum stress with yield stress, then for a given shell 
thickness, a Least weight solution will be possible, 

ff. The solution is found by solving the problem with frame 
area incremented a unit amount and with frame spacing incre- 
mented a unit amount. These two solutions yleld changes in 


the critical values of pressure and stress; that is, they 
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Give the partial differentials of pressure and stress 
With changes in frame area and frame spacing. It is 
then a simple matter to assume that linearization is 
possible and solve two equations in two unknowns for 
the number of unit changes to be made in frame area, A, 


and frame spacing, SP. 


oP oP 


— Jae = 2 

oA = oSP oe Foes Porit? 
OS os 

eee SSP co (Ss ela 7 Soi) 


solving for the two unknowns, dA and dSP, the original 
estimate of frame area and frame spacing can be changed. 
The problem is solved again with the new values of 

frame area and frame spacing and if the critical pressure 
is within a given percentage (the program specifies 4% 

on the safe side) and if the critical stress is also safe 
and within a set tolerance (these tolerances can be 
changed if a higher degree of accuracy is desired), the 
program will compute weight of compartment including shell 
and frames. (One word of warming. In order to ensure a 
satisfactory solution in all cases, it was necessary to 
make frame spacing continuous. Therefore, the number of 
frames, equal to compartment length over frame spacing 
may not be a whole number. It will be necessary to select 
the next higher number of frames to ensure a safe solution 


with a slightly higher weight than that computed.) 


= iO ee 





if, act the end of the first iteration, tne values of 
predicted collapse pressure and predicted stress are nov 
within tolerance, the program will continue to cycle 
through the serics of calculations improving the solution 
until it is within tolerance. The program then increments 
shell thickness to 95% of its predicted value and solves 
the problem again. The weight determined is compared with 
the weight determined with the original value of ‘t' and 
the "Least weight" is printed before the problem goes on 


to the next set of data. 


Seow 





LIST OF SYMBOLS 


some of the symbols in the following list refer to variables 


which are used to simplify the basic formulae. When such variables do 


not have explicit definitions, reference is made to a formula number, 


figure number, or page number, where the symbol is used in the notes 


by Professor Evans. 


cr 


CiiesoesecuLouo lenrame an ee 

Area of one frame plus shell under the web 
Web thickness of the ring frame 

Page 18 

Formula 31 

Outside diameter of hull 

Diameter of hull to mid-shell 

Young's Modulus 

Page 31 

Figure 8 

Moment of inertia of a T section 

Frame Spacing 

3.14 D ~ 2L (rounded off to the nearest whole number) 


2, 3, or 4 =~ no. of lobes for the three modes of 
failure by general instability 


Predicted collapse pressure 


Pressure at which yield occurs 


ow 





oS 


Lowest collapse pressure of the 6 pressures calculated 


Yield pressure - same as ze 
ressure at collapse deme. 
Radius of hull to mid-shell 


Highest stress of the three computed 


Frame spacing 

Shell thickness 

Page 34 

Page 34 

Page 18 

Page 9 

Slenderness ratio - Figure 4 


Yield stress 


Se 





a nput 


a. 


Input consists of the length and mean diameter of the 


cylinder under consideration, the characteristics of the 


material to be used in construction, including yieid 


strength, Young's modulus and density of the material, 


the collapse depth, allowable out-of-roundness of the 


hull, and the best estimate available of frame area and 


frame spacing. 


Inputs 
A 
CD 
CL 
D 
Delta 
Dens 
E 
Die 


OY 


Units 
ae 
feet 
inches 
inches 
inches 
1b/in? 
Weceire 


inches 


ie sells 


Definition 
Area of frame 
Collapse depth 
Compartment length 
Mean diameter of hull 
Allowable out-of-roundness 
Density of material 
Young's modulus 
Frame spacing 


Yield strength 


(Delta, allowable out-of-roundness of the hull, should 


normally be limited to one-half inch or half the thick- 


ness of the shell. 


For heavy shell plating this may be 


increased, but for "shell plauina Wess sudaneone—iice 


thick Delta will be altered by the program to one-half 


the thickness of the shell plating. ) 


~ kh - 








ibe 


All input for a single problem is placed on one card 


in the following form: 


The numbers indicate the card space number at which each item 


should start. Sample input follows: 


(0. ==--60000.--=-30000000.= 360. = = 00u= 5 4Oe a= eo Omar = 


Tae sample input indicates the temlewiae: 


Frame area TOMeQ= tire 

Yield Strength 60,000 pss me 
Young's Modulus 50s 10° Des eile 
Mean diameter of hull 360 inches 
Collapse depth 4,000 feet 
Compartment Length 540 inches 
Density of material -2875 lb. per in.> 
Frame spacing SOeime me. 


Allowable out-of-roundness 0.50 inches 
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An input deck ready for the computer will consist of the 


following cards: 


Meaning 


Identity 
Execute 
Binary 


Comment 


space 


Wumber 


7 
M3999-3504,FMS,DEBUG, 1, 1, 500, 75 Peterson 


XEQ 
BINARY (or FORTRAN ) 
Peterson (this card included only if a 


fortran deck is used) 


First binary or fortran card of program 


hast card jo: proeran 
END 
DATA 


First data card 


eee 






10 


20 


40 
30 


ED 


60 


64 


oe) 


67 


MINIMUM HULL WEIGHT PROGRAM 


DIMENSION PCAM(6)sPCW(6) sPCL(6),PCB11 6) sPCs2(ojePCe ws. 


1}»SMK1(6)s5MK2(6)5SMK3(6)s5P(10)35S04) swt 2G) er eee 
Zo Sis oul Oy) 


READ 1340 sAsSYsEsDsCDeCL»sDENSsSP»D=LTA 
PCD=CD*64./144. 

PRINT 1010,SY 

PRINT 10403CD 

PRINT 1090,PCD 

WSCD364. 4D 06S 7 Gl Cena Gy 

le tt =04.2 5.20 2 0 5.30 

T= One 5 

IF (T-120)40550950 

DELI A=172. 

CONTINUE 

SY¥=0s98%S Y 

CO2lw O02*CD 

PRINT 10205E 

PRINT 10303) 

PRINT 1050,CL 

PRINT 1060 ,DENS 

PRINT 1080 ,DELTA 

Pee woe 

L R= 

Lad 

Oe i 

LP=i 

NP=1 

N@=1 

PRELIMINARY CALCULATIONS 

FN=CL/SP+1l. 

IF (A) 64564565 

PRINT 1390,A 

GO TO 10 

IF (SP)66966567 

PRINT 140%, SP 

GO TO 10 

Qa i 

TW=SQRTF(A/19218) 

THETA=18.2*SP/ (D*¥SQRTF(100.*T/D)) 
SINH1=(2e71828**THETA=-1-e/2¢71828*¥*THETA)/2-0 
COSH1=(2e7T1B28*¥¥THETA+le/2Ze7T1E28%XTHETA! /220 
GK=(SINHI-SINF (THETA) )/(SINHI+SINFE( THETA) ) 
GN=(COSHI-COSF (THETA) 7 (oe ae a 
ire Het A 2, 
SINH2Z=(2-e71828**¥THA2Z—le/2e71828%*¥THA2Z)/2.6 
COSHZ=( 2.71828 #*THAZ+1e/2-71828**THA2)/2.6 
GH=-(36*SINH2*COSF (THA2)+COSH2*SINF(THA2))/(SINH1L+SINF 


L( THETA) ) 


Q==-(SINH2*COSF (THA2 )+COSH2*SINF(THA2Z))/(SINH1+SINF (THE 


ieee 


GT==$(SINA2*COSF ( THA2) -—COSH2RS TNE BmAZ yy (ol Ni oe eee 


amg g :. 





eo 
80 


LETA)) 
BDA=((SP/D)*#*2/(7T/D) #3) #206 25+50R eto) ee 
EPSL=(0.45*SY*D)/( EXT) 

B=TWHT/S(AtFTWXT ) 
BETA=(1l1lse#GN*T#*2)/USORIF (100.4070) Ch.) ae 
BM=3.1416204102---Gim) 

IEM=XFIXF (BM+0.5) 

BM=FLOATF(I8M) 

F=1.4+2.*Q*(0.85-B)/(1e+BETA) 

X=1e+GH*(0.85—-B)/(1¢+BETA) 
Z=0654+3e62*G6T* (0285-8) /(1.60+BETA) 
EL=1.57*S°RTF(D¥T7Z2.) 

IF (EL-SP)805803;70 

EL=SP 

ERT IAL=102.2*TWH*t4s(3.*%E LT / PEL eT +419 16 wore ee eee 
LT/TW+1 043) ¥#2+0.79) 

AB=A+TWHT 

C=(43 .5*#TWH*34(THEL) (8 OX TWH eS¥T))/S(19e1LOEXTWHE2Z4EL¥T ) 
PCAM(L)=2e*®T#SY/(D¥(16+GH¥((0-.85-B)/(1.e+BETA) ))*¥( 1040.6 
1148%*( (8BDA#*2—EPSL) #*2))) 


NO=1 
P(NO)=PCAM(L) 
NO=NO+#] 


PCOW(L)=2.*7T*#SY/(D*SQRIF (F442=F 7/2040, Zone) 

ONO y= Pew (CL) 

NO=NO+1 

PCL(L)=2,0*%T*SY/(D*SQRTF (X**2—-X*2Z24+2Z%*Z ) ) 

P(NO)=PCC(L) 

NO=NO+1 

BN=2. ; 
PCB1L(L)=2e* EXTHBM##4G/ (D* (BNHH#24+(2.5%BM) H¥2—-1.2)*%(BN¥#24+6 
IM#*#2)#%2)4+8e%(BN¥*¥2—-1.)X# EX ERTIAL/S (D*¥*3*SP) 
P(NO)=PCB1I(L) 

NO=NO+] 

BN=3e 

PCB2(L)=2e* EX T#ABME#X4/ (D¥(BNK¥24(25%BM) *¥¥*¥2—-1.)*(BN¥*24+3 
IM¥#2)¥*#2)+8e%(BN¥*X¥2~-1.)#EXERTIAL/S (D¥%3#SP) 
P(NO)=PCB2(L) 

NO=NO4+]1 

BN=4e 

PCB3(L)=2e* EX THBME*X4/ (D¥(BN¥*¥24+(.5*BM) ¥¥ 2-1.) *( BNX*24+5 
IM¥¥2 )*¥#2)4+8.%*(BN¥¥2—-1.)*E*XERTIAL/S(D*#3#SP) 
P(NO)=PCB3(L) 

Bit=Z24 

SMK1(L) =(PCD¥(005%D+065%T)*TW/ (AB¥(1,04+BETA)) )¥(1-0+(0 
eo eB ELA)/B) 
2+E¥C*%( BN*¥* 2-120) *DELTA*PCD/ ((0.5*D)**¥2%*(PCB1(L)—PCD) ) 
NO=1 

S(NO)=SMKRI(L) 

NO=NO+] . 

BN= 3 
SMK2(L)=(PCN¥(0.65*D4+0.5%F)#TW/(AB¥(1,0+4+BETA)) )¥(12¢0+(0 
ieee oebeTA)/ B) 
2+E%C#¥ (BN¥*¥ 2-120) *DELTA*PCD/((0.5*D)**¥2*(PCB1(L)-PCD) ) 
S(NO)=SMK2(L) 








ee 


20 


100 
£20 
eZ 2 
2 > 
130 


140 


ro 0 
CO 
6 2 
£3 
164 


1S 5 
70 


90 


NO=NO+] 

SN=4, 
SMK3(L)=(PCD*¥(0.5*D4+0.5%1)*TW/(AB*¥(1.0+BETA)) )*¥(1-0+(0 
1.85*BETA)/B) 
2+5%#C#(BNE*¥2—] 2.0) #DELTA*PCD/((0.5*D)**2%(PCBl1(L)=—PCD) ) 
S(NO)=SMK3(L) 

COMPUTE VOLUME OF MATERIAL IN ONE 7OMPARTMENT LENGTH 
V=0.7854*CL*( (047/26) #*#2-(D-1/26 ) 4x2) 405 (654% Dee 
1**2-(D-T/2.- 
25 693% TW) FEZ) 4 TWHEN+O 0. T8544FNX6 4HTWH( (D-T/2.-52693%1W) * 
3%2-(D-T/2.-8.0*TW) ¥¥2) 

COMPUTE WEIGHT OF MATERIAL IN ONE COMPARTMENT LENGTr 
W(L)=V*DENS 

L=L+1 

GO TO (905165) 5LP 

DO 110 N=156 

PDIFS(N)=PCD-—P(N) 

CONTINUE 

PDIFCR=PDIFS(1) 

Dem 20 N=256 

IF (PDIFS(N)=-PDIFCR) 12051003100 

PDIFCR=PDIFS(N) 

NP=N 

CONTINUE 

IF (0.02*PCD=$ABSF(PDIFCR) 113051253125 

LQ=LQ+] 

DO eee N=] ,3 

SDIFS(N) =SY-S(N) 

CONTINUE 

SDIFCR=SDIFS(1) 

Diemeke Oi Nl=2..3 

IF (SDIFS(N)-SDIFCR) 15091603160 

SDIFCR=SDIFS(N) 

NQ=N 

CONTINUE 

IF (0.02*SY-ABSF(SDIFCR) )16451633163 

LO=LQ+1 

GO MO (ieee 6 45.200) seo 

aie. 

PCRIT=PCU-PDIFCR 

SCRIT=SY-SDIFCR 

GO TO (17031803190) +LO 

Rs 

A=0296*A 

Ee-2 

GO TO 60 

DELPE1=PCRIT—P(NP) 

DELS1=SCRIT-S(NQ) 

AT 

SPT=SP 

SP=1.04*SP 

Lig 

GO TO 60 

DELPE2=PCRIT—P(NP) 

DELS2=SCRIT-S(NQ) 
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SP=o? | 
DOT=DELS2-DELPE2Z*DEU Say panei 
DOP=SDIECR=POUGGR*DEiS te aia 
PRY=DOP?/ Dom 
PRX=(FPOIFCR=PRY*DELPE2 ye eo 
A=(1e+0404*PRX)*#A 


200 


SP=(1 
GO TO 
Pie ill 
Ret | 
rs Neh 
iL INeL 
PRINT 
PRINT 
ea | 
PRINT 
PRINT 
eae 
PRINT 
eine) 
ae Nae 
PRINT 
PRINT 
a LNG 
PRINT 
PRINT 
aaa 
PRINT 
PRIN gt 
Pa ea 
ae Nel 
P ale te! 
PRIN | 
P RaniINet 
PRINT 
aie 
Pie 
PRINT 
PRINT 
Peri valt 
PRINT 
SN) 
PRINT 
Bia 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
See 
PRINT 
PRINT 
PRINT 
ele 
PRINT 


o-~eO4*PRY)#SP 


ts 

1000,A 
LOTOsFN 

1140 sPCAM(L=1) 
1150sPCW(L=1) 
LloOs»PCL(L=1 ) 
ily Ome esi =a} 
1200 ,PCB2(L-1) 
P21 Os PCS 4 (imi) 
12209SMK1(L-1) 
T2305 SMK2Z (lei) 
1240 »SMK3(L=1) 
1100%SP 

O51 

1120 5,7TW 

1180 5V 

1190 sW(L~-1) 
450sTHETA 
450sSINH1 

io Ome oi, 

450 GK 

450 GN 

450 5THA2 

MoO so) NZ 
BOO ee shi 2 
450%sGH 

4509Q 

4509G6T 
450>BDA 

Bo Og eRe 
45096 
450,sBETA 

450 3sBM 

450 99F 

4509X 

LOOm Z 

oO see 

BoOseERN LAL 
450sAB 

£505 C 

11305 PCR 
Peo Os oC Wa 

1 305 PD leer 


ai 30%5 
11305 
1130s 
eles U5 
Vi30% 


SOlrGr 
DELPEL 
Bieber 
Deo 
DE RSZ 
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220 
250 
240 


20 
260 


270 
450 
1000 
HO 10 
1020 
HOS 0 


1040 
eo > 0 


iO 7 0 


1060 


L080 
1090 
1100 
£110 
BL 20 
3 0 
1140 
ri) 0 
peL 60 
iO 


iL 30 


PRINT. 1.3042 poR 

PRINT 113050 ee 

PRINT 2.1305. ky 

er TN PS Oa eR 

GO 70 (2105220 s250)0 ae 

WLEAST=W(L=-1) 

PRINT 13705 WLEAST 

LR=2 

TR=T 

T=0-95%T 

GOs a> 

TF (WLEAST=W{( b=199 2450s 2405226 

WLEAST=W(L—1) 

PRINT. US7O,WeEAST 

T=1.04%TR 

LR=3 

GO Tees 5 

TF OWREA STIS CE Lyi 2 7082 7Os2o 

WLEAST=W(L=1) 

PRINT 21370, WLEAST 

T=1.04%T 

Go), Or 25> 

PRINT J1370,WLEAST 

GOm Oy VO 

mOmMAtt 0. 5) 

FORMAT( 1LH1LslLOHFRAME AREA219Xs1lH=sF10.233Xs7HSQe INe) 

FORMAT(CLIHILs14HYIELD STRENGTHs15XslH=5F100e2s 3X s6HPeSele) 

FORMAT(IXs14HYOUNGS MODULUS 5s12XslH=9F13e233Xs6HPeSele) 

FORMAT(1Xs21HMEAN DIAMETER OF HULL>s8X31lH=59F10e253X3s6H! 
NCHES) 

FORMAT(IXs14HCOLLAPSE DEPTHs15X3s1H325F10.293X54HFEET ) 

FORMAT(I1Xs16HBULKHEAD SPACINGs13X slH=sF1l00e2s3Xs6HINCHE 
5) 

FORMAT(1Xs30HNO. FRAMES BETWEEN BULKHEADS =3F10.2293X96 
HeRAMES) 

FORMAT(IXslOHDENSITY OF MATERIAL s1OX,1lH=sFl1lOe4s3Xs16HL 
Se reer CU 


Gr) 

FORMAT (1Xs30HALLOWABLE OUT-OF-ROUNDNESS = 9F100293X96 
HINCHES) 

PORMAT(IXs30HPRESSURE AT COLLAPSE DEPIF =9F10e293X%96 
iP eerie) 


FORMAT(1Xs13HFRAME SPACINGs:16Xs1lH=3sF100e233Xs6HINCHES ) 

FORMAT(1Xs15HSHELL THICKNESS314X91lH=5F10e23 3X s6HINCHES ) 

FORMAT(1Xe13HWEB THICKNESS 316X31lH=9F100233Xs6HINCHES ) 

FORMAT(F20-5) 

FORMAT (1HOs30HCOLLAPSE7PRESSe FROM 92AMOD. =3F100293Xs 
GHP Sele) 

FORMAT(1X»30HCOLLAPSE PRESS. FROM WENK =3F10.233X96 
EB Pis'S erties) 

FORMAT(1Xs30HCOLLAPSE PRESS FROM LUNCHICK 
ED’ Gerles) 

FORMAT(1X»s30HCOLL PRESS FR BRYANTs 2 LOBES=9F10e293X96 
Aiea Sees 

FORMAT(1Xs30HVOLUME OF ONE COMPT MATERIAL =sF100e233Xosl 


sF1l062s 3x86 
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ZAGUB I] Gather 


ioe on) 

1190 FORMAT(1Xs21HWEIGHT OF COMPARTMENT s8Xsln=sr 20.25. ocn. 
OUNDS) 

1200 FORMAT(1Xs30HCOLL PRESS FR BRYANIS) 3S LOse>=5F Ogee ae 
HP eSele) 

1210 FORMAT(1X»30HCOLL PRESS FR BRYANT S 4 COSES=— 290 ane 
HPeeSele) 

1226 FORMAT(1Xs30HMAX STR FR BUSHIPSs 2 LOBES =sFlO e255 
Tie Sec! on) 

1230 FORMAT(1X»,30HMAX STR FR BUSHIPSs 3 LOBES =J 16.2. 
HP eSele) 

1240 FORMAT(1Xs30HMAX STR FR BUSHIPSs 4&4 LOBES =59F10e233X96 
HP eSele) 


mo4+0O FPORMATUTFGO*2s2F104033F6.035F Gees ee esr. 
13760 FORMAT (IX s30HCOMPARTMENT OF LEAST WEIGHT =9F10.233Xs6 


HPOUNDS) 

1390 FORMAT(1Xs30HACTION TERMINATED, NEGATIVE A=sF1l0e233X97 
Sere | INcew) 

1400 FORMAT(1Xs30HACTION TERMINATED, NEGATe SP =sF10e233X26 
PpuINGre oo 


END 


NOTE 


FORMAT STATEMENTS REQUIRE A CARD FORM 72 CHARACTERS 
IN WIDTH. IN ORDER TO FACILITATE A NARROWER PRINTOUT 
mROSE FORMAT STATEMENTS LONGER THAN ©6O CHARACTERS eke 
BROKEN AND THE EXTRA CHARACTERS PRINTED ON THE FOLLOWING 
LINE WITH AN EXAGERATED INDENTATION. A USER COPYING THE 
PROGRAM SHOULD TYPEs FOR EXAMPLE, LINE 1400 AND THE NEXT 
INDENTED LINE ON THE SAME CARD. 
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BDA 


BETA 
BN 


CD 


COSHL 


COSH2 
D 
DELPEL 
DELPE2 
DELSL 
DELS2 
DELTA 
DENS: 


DOT, DOP 


COMPUTER PROGRAM SYMBOLS 


Cross-sectional area of a frame 


Cross-sectional area of a frame plus the area of 
Shell under the web 


pame as A 

Page 18 

Lambda - Figure 4 

Page 18 

m (see formula symbols) 
n (see formula symbols ) 
Formula 31 

Collapse depth 
Compartment Length 
COSH © 

COSH 6/2 

Hull diameter 


P - Pp (after A is incremented) 


eri ert 
Bes eae (after SP is incremented) 
Sorit ~ eer (after A is incremented) 
Seg > Sees (after A is incremented) 


Allowable out-of-roundness of the hull 
Density of construction material 


Intermediate values in the solution of two 
differential equations 


Young's modulus 
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= 


EPSL 


ERTIAL 


rf 


PN 


GN 


GT 


PCAM 


PCBL 


ECre 
PCB3 
PCD 
PCL 
PCRIT 
PCW 
PDIEFCR 
PDIFS 


PRX 


ERY 


Kffective length of shell plating 

e, page 9 

Moment of inertia of a T section 

Page 31 

No. of frames in one compartment length 
Figure 8 

Figure 8 

Figure 8 

Page 34 


Predicted collapse pressure = PCAM, PCW, PCL, PCBL, PCBe, 
PCBS 


Predicted collapse pressure from modified formula 92A 


Predicted collapse pressure from Bryant's (Kendrick) 
formula (lst failure mode) 


(2nd failure mode) 

(3rd failure mode) 

Pressure at collapse depth 

Predicted collapse pressure from Lunchick 
Lowest collapse pressure of the 6 predicted 
Predicted collapse pressure from Wenk 
Highest of the 6 PDIFS values 

= PCD - P (six of these calculated ) 


No. of increments by which SP must be changed in 
reaching an optimum solution 


No. of increments by which A must be changed in reaching 
an optimum solution 


Page 31 


- 54 - 





SCRIT 
SDIFCR 
SDIFS 
SINHL 
SINH2 
OMK1L 


SMK2 


OP 


SPT 


oY 


HAC 


THETA 


WLEAST 


Predicted stress = SM, sMhz mes 
Highest stress of the 3 predicted 

Lowest of the 3 SDIFS values 

SY - S (three calculated) 

SINH © 

SINH 9/2 

Predicted stress (2 lobe mode of failure) 
Predicted stress (3 lobe mode of failure) 
Predicted stress (4 lobe mode of failure 
Frame spacing 

pame as SP 

Yield stress 

Shell Thickness 

@/2 

© (Figure 8) 

Web thickness 


Volume of material in shell and frames over one 
compartment length 


Weight of material in shell and frames over one 
compartment Length 


Weight of the lightest compartment of those considered 
Page 34 


Page 34 
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Le Output 


ae 


Given the following inputs for a section of pressure 


mules 


1.) Diameter 

2.) Length between transverse bulkheads 
3.) Collapse depth 

4.) Construction material 


5.) Tolerance to which the construction 
material will be fabricated 


the computer program produces recommended values for: 
1.) Shell thickness 
2.) Frame size, and 


3.) Frame spacing. 


Wit the above information, it is possible to compute 
displacement of the cylinder and weight of material used 
in construction of the hull. The computer program does, 
in fact, calculate weight of hull since it was decided 
that “least weight" would be the criterion for optimum 
design. 

But the process of design involves more than one answer 
to one problem. Normally, a designer will study several 
vehicles with various characteristics and he will attempt 
to gather as much information as possible about each one. 
It is not enough to know that a computer program has 


recommended certain values for a certain vehicle. Each 
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Ot ene recommended values must be checked for accuvacwe 
logic, and applicability to the »lannea use end the 
construction mevhods currently available. There is 
always the possibility that a computer program will 
recommend as optimum a set of values impossible to adhere 
weetieacuual Construction. 

To aid a designer in comparing vehicles and to aid in 
Qralyzing tae computer program itsely, ox possibly te 
help in studying the design parameters of submarine con- 
struction, the computer program listed in this section 
prints out a large number of values computed in the 
preliminary calculations as well as the final answers 
to the main problem. Once a user has gained familiarity 
with the program he will probably eliminate much of the 
printed output by withdrawing from the input deck those 
PRINI statements referring to the unnecessary values. 

He may also add some cards or some additional calculations 
which will help him in his design or analysis. 
Referring to the sample output (page 60), the following 


values are printed out after each solution av a. given 


thickness. 
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Computer Output 


ike 


3-8. 


9-11. 


ee 
13. 
ad 
15. 
MO 
7-38. 


39-906 


Dilue 


Page One 


Yield strength of the material 


Collapse depth 

Pressuce au. Collapse depth 

Young's modulus 

Hull diameter 

Compartment Length 

Density of the material 

Allowable out-of-roundness of the hull 
Page Two 

Frame Area 

Number of frames 

Critical pressure predicted by formula 92a- 
modified, Wenk, Lunchick, and three modes 

by Bryant 

Three values of yleld strength predicted by the 
formula in section l-c-6 of this appendix. 

Frame spacing 

shell thickness 

Web thickness 

Volume of material in the shell and frames 

Weight of material in the shell and frames 
Preliminary calculation values beginning with "THETA" 
and ending with "C" as Listed on the fourth page of the 
computer program. 

Values used in the solution of the two differential 
equations when determining the required changes. 
"Least weight" of those compartments considered. 


sone 





d. Since thickness is first decreased to 95% of the 
predicted value and then increased to 104% of that 
predicted by the hoop stress formula, the values on 
page two will be printed at least three times. IT, 
as may happen in some cases, the solution with 
thickness increased produces the "least weight" 
Solution, thickness will continue to be increased by 
increments of 4% until a "least weight" solution is 
reached and passed. (Each of these incremented 
Solutions will also require one page of printed 
Output.) his Series of events wouldyoceu: if sib), 
increasing thickness, frame size and the number of 
frames could be decreased enough to achieve an 
overall reduction in weight. After a certain number 
of thickness increases, however, weight will begin 
to increase at which time the program will terminate. 
Experience indicates, however, that the "least 
weight" solution will occur when "I" is reduced to 
95% of the thickness predicted by the hoop stress 
Homma, and if the hoop stress prediction 1s igueveds 
"least weight" will oceur at about 70% of the 


predicted thickness. 
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